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Monte Carlo simulations of the two-dimensional XY model are performed in 
a square geometry with free and mixed fixed-free boundary conditions. Using a 
Schwarz-Christoffel conformal mapping, we deduce the exponent 77 of the order pa- 
rameter correlation function and its surface equivalent 77|| at the Kosterlitz-Thouless 
transition temperature. The well known value ^(Tkt) = 1/4 is easily recovered 
even with systems of relatively small sizes, since the shape effects are encoded in the 
conformal mapping. The exponent associated to the surface correlations is similarly 
obtained 771 (T K t) — 0.54. 

1 Introduction 



The two-dimensional classical XY 
model exhibits a non standard behaviour. 
Long-range order of spins with continuous 
symmetries is indeed forbidden according 
to the Mermin- Wagner theorem |ffl], and 
thus there is no spontaneous magnetisa- 
tion at finite temperature. The transition 
is governed by unbinding of topological de- 
fects. The low-temperature regime was in- 
vestigated by Berezinskii |2| and the mech- 
anism of unbinding of vortices was studied 
by Kosterlitz and Thouless || using renor- 
malization group methods. For reviews, see 
e.g. Refs. |,|. 

The notations are specified below: we 
consider a 2d XY-model on a square lattice 



with two-components spin variables (fi w = 
(cos6 w ,sin9 w ) located at the sites w of & 
lattice A of linear extent L. The spins in- 
teract through the usual nearest-neighbour 
ferromagnetic interaction 



H 



(1) 



where K = J/ksT and £1 is a unit vector 
in the fi— direction. 

At low temperature, the system is par- 
tially ordered, apart from the existence of 
vortices, which appear in pairs in increas- 
ing number with increasing temperature. 
In the low-temperature limit, the effect of 
vortices can be neglected and the behaviour 
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is governed by spin- wave excitations, ob- 
tained after expanding the cosine. This 
harmonic approximation is justified, pro- 
vided that the spin disorientation remains 
small, i.e. at sufficiently low temperature. 
Within this approximation, the quadratic 
energy in the Boltzmann factor leads to a 
Gaussian equilibrium distribution and the 
two-point correlation function becomes H 

' < -1/2** (2) 



\W\ — U>2\ 



hence 

r)(T) = k B T/27rJ. (3) 

When the temperature increases, the influ- 
ence of vortices becomes more prominent, 
producing a deviation from the linear spin- 
wave contribution in equation (Q), but the 
order parameter correlation function still 
decays algebraically with an exponent r](T) 
which depends on the temperature. The 
existence of such a scale-invariant power- 
law decay of the correlation function is 
the signature of a continuous line of fixed 
points at low temperatures. In this paper 
we are interested in the end point of this 
critical line where the topological transition 
takes place: at the transition temperature 
Tkt (usually called Berezinskii-Kosterlitz- 
Thouless critical temperature), the pairs 
are broken and the system becomes com- 
pletely disordered. This very peculiar topo- 
logical transition is characterised by essen- 
tial singularities when approaching the crit- 
ical point from the high temperature phase, 
and at Tkt the correlation function expo- 
nent takes the value ^(Tkt) = 1/4 ||. 



In order to recover this result, many nu- 
merical studies were performed in the last 
decade using Monte Carlo simulations (see 
e.g. Refs. 



1C, 11, 12 ), but the anal- 



ysis was made difficult by the existence of 
logarithmic corrections, e.g. 



hi2r"(ln M)^ 2 , (4) 



at Tkt ||- Due to these logarithmic cor- 
rections which make the fits quite difficult, 
the value of r/ at Tkt was a bit controversial 
as shown in table 1 of reference [jT^] . The 
resort to large-scale simulations, up to sys- 
tems as large as L = 1200, was then needed 
in order to confirm this picture [11]. 

Recently, we proposed a rather differ- 
ent approach |l3|] which is easily imple- 
mented. We use the covariance law of 
correlation functions under the mapping 
of a two-dimensional system confined in- 
side a square onto the half-infinite plane. 
The scaling dimensions are then obtained 
through a simple power-law fit when the 
correlators are expressed as functions of 
conveniently rescaled variables. The shape 
effects are explicitly encoded in the con- 
formal mapping. This is the crucial point, 
since then even very small systems lead to 
promising results. We may then determine 
accurately the bulk correlation function ex- 
ponent r] at Tkt but also, by choosing con- 
venient boundary conditions, the surface 
critical exponent which describes the decay 
of the correlation function parallel to a free 
surface, m. 



2 Critical system confined inside a square with open 
boundaries 

For a scale- invariant system at its crit- teractions), conformal invariance enable to 
ical point (which also exhibits the proper- include explicitly the shape dependence in 
ties of isotropic scaling and short-range in- the functional expression of the correlators 
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through the conformal covariance transfor- 
mation under a mapping w(z): 

{$ wi 4w2) = \w\zi)\~ Xa \w\^)\~ x °{$zx4>zz)- 

(5) 

The functional expression of the correla- 
tion function inside the square geometry w 
simply follows the mapping w(z) which re- 
alizes the conformal transformation of the 
half-plane z = x + iy(0<y< oo) in- 
side a square w = u + iv of size L x L 
(-L/2 < u < L/2, < v < L) with free 
boundary conditions along the four edges. 
This is realized by a Schwarz-Christoffel 
transformation [^] 

w(z) = ^F(z,k), 2 = snf^-^.(6) 

Here, F(z, k) is the elliptic integral of the 
first kind, sn (2Kw/L) the Jacobian ellip- 
tic sine, K = K(k) the complete elliptic 
integral of the first kind, and the modu- 
lus k depends on the aspect ratio of A and 
is here solution of K(k)/K(\Jl — k 2 ) = |. 
Using the mapping @, one obtains the lo- 
cal rescaling factor in equation (|5|), w'{z) = 



2^[(1 - z 2 )(l - k 2 z 2 )}- 1 / 2 . One simply has 
now to include in equation (||) the expected 
behaviour of (<j) zi ■ (p Z2 ) in the half-infinite 
geometry in order to get the functional ex- 
pression of the corresponding quantity in 
the confined system. For example in the 
semi-infinite geometry z = x + iy (the free 
surface being defined by the x axis), the 
two-point correlation function is fixed up 
to an unknown scaling function: fixing one 
point z\ close to the free surface {z\ = i) of 
the half-infinite plane, and leaving the sec- 
ond point Z2 explore the rest of the lattice, 
the following behaviour is expected: 

(4 • K) ioo ~ (ui ~ y2)~ Xa ^), (7) 

with T) = 2x a . The dependence on u = 
i yiV2 m of the universal scaling function 

|zi— Z2| 

i\) is constrained by the special conformal 
transformation and its asymptotic be- 
haviour is implied by scaling e.g. ip(uj) ~ 
oj x t when y2 3> 1, with x\ = ±7/m the sur- 
face scaling dimension. Keeping w\ fixed 
inside the square, the two-point correlation 
function becomes 



k(w) 



»w)f 

- Im 



[k(w)\ Xa ip(uj) 
2Kw~ 



sn- 



L 



1 



2 2Kw\ ( 



where / specifies that the boundary condi- 
tions (BCs) are chosen free along the edges 
of the square. This expression is correct 
up to a constant amplitude determined by 
k(w%) which is kept fixed, but the function 
ip(uj) is still varying with the location of the 
second point, w. This makes this expres- 
sion not so useful in practice (see e.g. ref- 
erence [17 1 for an application). In order to 



cancel the role of the unknown scaling func- 
tion, it is more convenient to work with a 



kV 2 -fL) 



-1/2 



(8) 



density profile m(w) in the presence of sym- 
metry breaking surface fields hg\ on some 
parts of the boundary <9A of the lattice A. 
In the case of fixed-free BCs denoted by +/ 
(the spins located on the x > half axis are 
kept fixed while those of the x < are free) , 
Burkhardt and Xue have shown in the case 
of Ising and Potts models that in the half- 
plane the order parameter profile obeys the 



following expression [18]: 



m +f (z) =y- r il 2 {cos \ef\\l 2 . (9) 
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In the square geometry, the fixed-free BCs u > 0, and the corresponding profile is ex- 
correspond to keeping the spins fixed for pected to obey the following ansatz: 



L \ ~ XrT i 
m+ f (w) ~ ( — J [k(w)]- x ° [h(w)] x ° , 



M») = 4= U + Re 



2Kw 



sn- 



L 



2Kw 



sn- 



L 



1/2 



(10) 



where k(w) was defined in equation @. 

In this paper we will essentially apply 
equation (|l0|) in order to determine the 
critical exponents rj and 7711 at the KT point. 



In Ref . |13[ , the bulk exponent was deduced 
from a similar analysis of the order parame- 
ter profiles in a square-shaped system with 
fully fixed BCs. 



3 Monte Carlo simulations 



The application of the simple power- 



law in equation (10) requires a relatively 
precise numerical determination of the or- 
der parameter profile of the 2d XY model 
confined inside a square with mixed fixed- 
free boundary conditions. Fixing the spins 
along some part dA of the boundary, e.g. 
4> w = (1,0), \/w £ <9A, plays the role of 
ordering surface fields hg^ w y The resort 
to cluster update algorithms (Wolff algo- 
rithm here) is necessary in order to pre- 
vent the critical slowing down: the cen- 
tral idea is the identification of clusters of 
sites using bond variables connected to the 
spin configuration. The spins of the clus- 
ters are then updated and the algorithm 
is particularly efficient if the percolation 
threshold of the bond process occurs at the 
transition temperature of the spin model. 
This ensures the updating of clusters of 
all sizes in a single MC sweep. For the 
XY model, Ising variables defined in the 
Wolff algorithm by the sign of the projec- 
tion of the spin variables along some ran- 
dom direction. The bonds are then intro- 
duced through the Kasteleyn-Fortuin ran- 
dom graph representation. The percolation 



threshold for these bonds coincides with 
the Kosterlitz-Thouless point [20], which 



guarantees the efficiency of the Wolff clus- 
ter updating scheme || at Tkt- When 
one uses particular boundary conditions, 
with fixed spins along some surface as re- 
quired here, the Wolff algorithm should be- 
come less efficient, since close to critical- 
ity the unique cluster will often reach the 
boundary and no update is made in this 
case. To prevent this, we use the symme- 
try properties of the Hamiltonian (|l|): even 
when the cluster reaches the fixed bound- 
aries, it is updated, and the order param- 
eter profile is then measured with respect 
to the new direction of the boundary spins, 
m(w) = (<fi w ■ (/>dA)sq. • The new configura- 
tion reached has the right statistical weight. 

Using this procedure, we discard 10 6 
sweeps for thermalization, and the mea- 
surements are performed on 10 6 production 
sweeps. The simulations at Tkt on a typ- 
ical system of size 100 x 100 takes a few 
hours on a standard PC with 733 MHz pro- 
cessor, to get a smooth profile as shown in 
figure 1. 
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Figure 1: MC simulations of the 2d XY model 
inside a square 48 x 48 spins with mixed fixed- 
free boundary conditions at the KT transition 
temperature IzbTkt/J — 0.893 (average over 
10 6 MCS/spin after cancellation of 10 6 for ther- 
malization). 



4 Scaling of the order parameter profile m + f(w) 



In Ref. [13], we have shown that the 
logarithmic corrections at the KT point 
are negligible for the order parameter pro- 
file with fixed BCs. Assuming that there 
will be no substantial difference for fixed- 
free BCs, we perform a fit of lnm+/(u;) 
against the two-dimensional linear expres- 
sion, const — Xfj In n(w) + x\ In fj,(w), as ex- 
pected from equation (fu]). The value ob- 
tained for the bulk correlation function ex- 
ponent rj will be a test for the absence of 
logarithmic corrections. The simulations 
are performed at the value ksT^/J ~ 
0.893 found in the literature, e.g. in 
Ref. ||. The fit leads to accurate deter- 
minations of both exponents for different 
system sizes ranging between L = 20 and 
L = 200. Although it would easily be pos- 
sible to produce simulations for larger sys- 
tems, there is no need here, since the results 
are already quite stable as shown in table 1 
and in figure 2. 



Table 1: 

Bulk and surface scaling dimensions of the or- 
der parameter at the Kosterlitz-Thouless tran- 
sition. The values that we quote are those 
which correspond to a fit over all the lattice 
points (except those exactly at the boundaries) 
ignoring possible lattice effects. 



L 




< 


xVd.o.f 




20 


0.149(33) 


0.290(15) 


0.691 x 10" 


-3 


24 


0.147(34) 


0.287(15) 


0.719 x 10" 


-3 


28 


0.134(22) 


0.280(12) 


1.210 x 10" 


-3 


3G 


0.133(35) 


0.278(12) 


1.010 x 10" 


-3 


48 


0.133(35) 


0.279(10) 


0.743 x 10" 


-3 


64 


0.131(35) 


0.278(10) 


0.664 x 10" 


-3 


100 


0.126(17) 


0.275(14) 


0.544 x 10" 


-3 


200 


0.125(23) 


0.271(11) 


0.392 x 10" 


-3 



This is fully coherent with the result of 
Kosterlitz and Thouless RG analysis r\ = 
2x a = 1/4. We can thus also consider as 
reliable the value of the surface exponent 

7711 = 2x1 — 



0.4 r 



0.3- 



i I * 



0.2- 



,u H 



0, 



0.01 0.02 



0.03 
ML 



0.04 0.05 



Figure 2: Size dependence of the bulk and 
surface scaling dimensions at the KT transition 
plotted as a function of the inverse linear size 
of the lattice. 
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5 Scaling of the correlation function (<j) Wl ■ <p w )f with free 
boundary conditions 



Equation @ provides a test in order 
to control the quality of the surface expo- 
nent obtained above. Indeed, if one plots 
{4> Wl • 4>w)sq. x [/«(io)] a:o ' vs cj, one should 
obtain the universal scaling function 
whose asymptotic large distance behaviour 
should lead to a consistent surface expo- 
nent through the expected power-law be- 
haviour oj x ° . 




6 Conclusion 

In this paper we have performed stan- 
dard MC simulations of the 2d XY model 
at the Kosterlitz-Thouless topological tran- 
sition. Combined to the use of a conformal 
mapping, they provide a simple numerical 
determination of bulk and surface critical 
exponents associated to the algebraic de- 
cay of the correlation function. The well 
known result of Kosterlitz and Thouless for 
the bulk exponent is recovered and as far as 
the author knows this is the first determi- 
nation of the surface exponent 7711 ~ 0.54. 
On the contrary to the bulk case, 7711 is dif- 
ferent from the Ising value (which is equal 

to ^ = 1) m. 



Figure 3: Universal scaling function of the 
surface-volume correlation function for L = 
100. The power-law behaviour tp(uj) ~ u x ° 
shown in dotted line agrees with a value x\ ~ 
0.277. 

The scaling function i/j(uj) is shown in 
figure 3. It exhibits a clear power-law be- 
haviour consistent with the previous value 
of the surface scaling dimension. 

We mention here that the data were 
obtained for both quantities, {4> Wl • 4> w ) f 
and m+f(w), with completely independent 
runs, since the boundary conditions are dif- 
ferent in equations (B) and ( |i"0| ) and the 
quantities themselves are different (2-point 
or 1-point correlators). The data collapse 
is quite good and supports the validity of 
equation (^). 
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